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DESCRIPTION 

NOx DISCHARGE QUANTITY ESTIMATION METHOD 
FOR INTERNAL COMBUSTION ENGINE 

TECHNICAL FIELD 

The present invention relates to an NOx discharge quantity 
estimation method which is applied to an internal combustion engine 
equipped with an EGR apparatus for circulating to an Intake passage a 
portion of exhaust gas flowing through an exhaust passage, and which is 
adapted to estimate the quantity of NOx contained in the exhaust gas 
discharged from the exhaust passage to the outside. 

BACKGROUND ART 

In an internal combustion engine such as a spark-Ignition engine or 
a diesei engine, the quantity of NOx contained in exhaust gas discharged 
from an exhaust passage to the outside (hereinafter may be referred to as 
"NOx discharge quantity") must be reduced. An effective way of reducing 
the NOx discharge quantity is lowering the highest flame temperature 
(highest combustion temperature) through, for example, increasing the 
quantity of EGR gas circulated by means of an EGR apparatus, or delaying 
fuel injection timing. 

However, when the quantity of EGR gas is increased in order to 
reduce the NOx discharge quantity, in the case of a diesei engine, the 
generation quantity of particulate matter (PM) increases. When fuel 
injection timing is delayed in order to reduce the NOx discharge quantity. 
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fuel efficiency deteriorates. 

Accordingly, in order to minimize tlie NOx discharge quantity in 
consideration of suppression of an increase in tiie discliarge quantity of 
particulate matter (PM) and suppression of deterioration in fuel efficiency, 
tlie NOx disctiarge quantity is desirably controlled to a predetermined target 
value corresponding to the operating conditions of the engine. In order to 
accurately control the NOx discharge quantity to a predetermined target 
value, the NOx discharge quantity must be accurately estimated. 

For such accurate estimation, a control apparatus for an internal 
combustion engine disclosed in Japanese Patent Application Laid-Open 
{kokaii No. 2002-371893 detects combustion pressure and intal<e-gas 
oxygen concentration by use of a cylinder pressure sensor and an 
intal<e-gas oxygen concentration sensor, and estimates the quantity of NOx 
generated as a result of combustion (hereinafter referred to as 
"combustion-generated NOx quantity") on the basis of combustion 
temperature and gas mixture concentration calculated on the basis of the 
combustion pressure and the intal<e-gas oxygen concentration, wherein the 
estimation is performed by use of the extended Zeidovich mechanism, which 
is a typical known combustion model. Then, EGR gas quantity, fuel 
injection timing, or the lil<e is controlled so that the estimated 
combustion-generated NOx quantity coincides with a predetermined target 
combustion-generated NOx quantity, whereby the NOx discharge quantity is 
controlled to a predetermined target value. 

Incidentally, in an internal combustion engine equipped with an EGR 
apparatus, NOx contained in EGR gas is circulated into a combustion 
chamber via the EGR apparatus. In addition, the above-mentioned 
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combustion-generated NOx quantity is the quantity of NOx generated in tlie 
region wliich is a portion of tlie combustion ciiamber and in whicii 
combustion occurs (iiereinafter referred to as "combustion region"). 
Accordingiy, in tlie remaining portion of tiie combustion ciiamber 
(iiereinafter referred to as "non-combustion region"), the circulated NOx 
remains after combustion. Therefore, in order to accurateiy estimate the 
NOx discharge quantity, not only the combustion-generated NOx quantity but 
also the "quantity of NOx remaining in the non-combustion region" must be 
taken into consideration. 

However, in the above-mentioned conventional apparatus, the 
"quantity of NOx remaining in the non-combustion region" Is not taken into 
consideration at all. Accordingiy, when the conventional apparatus Is 
applied to an internal combustion engine equipped with an EGR apparatus, 
NOx discharge quantity cannot be accurately estimated, and thus, the NOx 
discharge quantity cannot be accurateiy controlled to a predetermined target 
value. 

DISCLOSURE OF THE INVENTION 

in view of the foregoing, an object of the present invention is to 
provide an NOx discharge quantity estimation method used with an internal 
combustion engine equipped with an EGR apparatus, which method can 
accurateiy estimate the quantity of NOx discharged from an exhaust 
passage to the outside, in consideration of quantity of NOx remaining in a 
non-combustion region. 

The NOx discharge quantity estimation method according to the 
present invention is applied to an internal combustion engine equipped with 
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an EGR apparatus, and adapted to estimate a quantity of NOx contained In 
exhaust gas discharged from an exhaust passage of the engine to the 
outside (i.e., the above-mentioned NOx discharge quantity) in consideration 
of a quantity of NOx generated in a combustion chamber as a result of 
combustion (i.e., the above-mentioned combustion-generated NOx quantity) 
and a quantity of NOx circulated Into the combustion chamber via the EGR 
apparatus. 

The quantity of NOx circulated into the combustion chamber via the 
EGR apparatus (partially) serves as the above-mentioned 'quantity of NOx 
remaining in the non-combustion region" after combustion. Accordingly, 
under the method of the present Invention, the NOx discharge quantity can 
be estimated in consideration of not only the combustion-generated NOx 
quantity, but also the "quantity of NOx remaining in the non-combustion 
region." Therefore, the NOx discharge quantity can be accurately 
estimated. 

A more specific method for estimating NOx discharge quantity 
according to the present invention, which is applied to an internal 
combustion engine equipped with an EGR apparatus, comprises the steps 
of estimating a combustion region, the combustion region being a region of 
the combustion chamber In which combustion occurs; estimating, as a 
combustion-generated NOx quantity, a quantity of NOx generated In the 
combustion region as a result of combustion, and an NOx quantity In a 
non-combustion region, the non-combustion region being the remaining 
region of the combustion chamber; and estimating, on the basis of the 
combustion-generated NOx quantity and the NOx quantity in the 
non-combustion region, a quantity of NOx contained in exhaust gas 
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discharged from the exhaust passage to the outside. 

Under this method, the NOx discharge quantity can be estimated in 
consideration of not only the combustion-generated NOx quantity; i.e., the 
quantity of NOx which is generated in the estimated combustion region as a 
result of combustion, but also the NOx quantity in the non-combustion region 
(after combustion); i.e., the above-mentioned "quantity of NOx remaining in 
the non-combustion region," Thus, the NOx discharge quantity can be 
accurately estimated. 

in this case, the NOx quantity in the non-combustion region to be 
estimated is preferably a non-combustion-reglon circulated NOx quantity 
which represents a quantity of a portion of NOx circulated into the 
combustion chamber via the EGR apparatus, the portion of the circulated 
NOx being present in the non-combustion region before combustion. 

Combustion within the combustion chamber occurs only in the 
combustion region, and does not occur in the non-combustion region. 
Therefore, of the quantity of NOx circulated into the combustion chamber via 
the EGR apparatus, 

the quantity of NOx present within the non-combustion region before 
combustion (before the start of combustion) (hereinafter referred to as 
"non-combustion-region circulated NOx quantity") can be considered to be 
conserved (held) as is after combustion. In other words, the 
non-combustion-region circulated NOx quantity can be used as the "quantity 
of NOx remaining in the non-combustion region." Accordingly, through 
employment of the above-scheme, the "quantity of NOx remaining in the 
non-combustion region" can be accurately estimated by use of a simple 
configuration (combustion model). As a result, the NOx discharge quantity 
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can be accurately estimated through simple calculation. 

Moreover, the combustion-generated NOx quantity Is preferably 
estimated in such a manner that when the combustion-region circulated NOx 
quantity is greater than the combustion-generated NOx quantity, the 
combustion-region circulated NOx quantity Is employed as the 
combustion-generated NOx quantity (which serves as the quantity of NOx 
remaining in the combustion region after combustion), wherein the 
combustion-region circulated NOx quantity represent a quantity of a portion 
of NOx circulated Into the combustion chamber via the EGR apparatus, the 
portion of the circulated NOx being present in the combustion region before 
combustion. 

Before start of combustion, NOx circulated into the combustion 
chamber via the EGR apparatus (i.e., the above-mentioned 
combustion-region circulated NOx quantity) is present in the combustion 
region, and after combustion, NOx generated as result of combustion (i.e., 
the above-mentioned combustion-generated NOx quantity) is added thereto. 
In general, nitrogen molecules N2 and NOx (In particular, carbon monoxide 
NO) are substances which can generate each other through reversible 
reaction. Therefore, the NOx quantity (NOx concentration) of a gas after a 
chemical reaction is determined, irrespective of NOx quantity (NOx 
concentration) before the chemical reaction, by means of chemical 
equilibrium, In accordance with the temperature of gas after the chemical 
reaction (in such a manner that the higher the gas temperature, the greater 
the NOx quantity (NOx concentration)). 

Moreover, in an ordinary operating state of the engine, in many 
cases the combustion-generated NOx quantity (quantity of NOx generated 
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within gas having an increased temperature as a result of combustion) is 
greater than the above-described combustion-region circulated NOx quantity. 
Accordingly, basically, the quantity of NOx generated as a result of 
combustion can be considered to be substantially equal to the 
combustion-generated NOx quantity (which serves as the quantity of NOx 
remaining in the combustion region after combustion). 

However, in a case where the NOx concentration of EGR gas is 
considerably high and fuel injection quantity is very small because of a very 
small accelerator opening, the combustion-region circulated NOx quantity 
increases, and the combustion-generated NOx quantity decreases because 
of a small Increase In gas temperature stemming from combustion. As a 
result, the combustion-region circulated NOx quantity becomes greater than 
the combustion-generated NOx quantity in some cases. In such a case, 
the combustion-region circulated NOx quantity, rather than the 
combustion-generated NOx quantity, is considered to be a more accurate 
representation of the combustion-generated NOx quantity (which serves as 
the "quantity of NOx remaining in the combustion region after combustion"). 

In other words, the combustion-generated NOx quantity (which 
serves as the quantity of NOx remaining in the combustion region after 
combustion) is preferably estimated by selecting the larger one of the 
combustion-region circulated NOx quantity and the combustion-generated 
NOx quantity (generated after combustion). Accordingly, through 
employment of the above scheme, the combustion-generated NOx quantity 
(which serves as the quantity of NOx remaining in the combustion region 
after combustion) can be accurately estimated by use of a simple 
configuration (combustion model). As a result, the NOx discharge quantity 
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can be accurately estimated through simple calculation. 

In the NOx discharge quantity estimation method according to the 
present invention, the NOx discharge quantity is preferably estimated by the 
steps of estimating an NOx concentration of exhaust gas on the basis of the 
combustion-generated NOx quantity and the NOx quantity in the 
non-combustion region; and estimating the NOx discharge quantity by 
multiplying the NOx concentration by a quantity of the exhaust gas 
discharged from the exhaust passage to the outside. 

The NOx concentration of exhaust gas discharged from the engine to 
the exhaust passage via the exhaust valve can be easily calculated as the 
ratio, to the total quantity of gas within the combustion chamber, of the "sum 
of the combustion-generated NOx quantity and the NOx quantity in the 
non-combustion region." Moreover, when the NOx concentration of 
exhaust gas flowing through the exhaust passage is assumed to be constant 
throughout the entire region of the exhaust passage, the NOx concentration 
of exhaust gas discharged from the engine to the exhaust passage via the 
exhaust valve becomes equal to the NOx concentration of the exhaust gas 
discharged to the outside from the exhaust passage. 

Accordingly, as described above, the NOx concentration of exhaust 
gas (at the time of discharge from the engine to the exhaust passage via the 
valve), is calculated on the basis of the combustion-generated NOx quantity 
and the NOx quantity in the non-combustion region, and the NOx 
concentration is multiplied by the quantity of the exhaust gas discharged 
from the exhaust passage to the outside. Thus, the NOx discharge quantity, 
the quantity of NOx discharged from the exhaust passage to the outside, can 
be accurately estimated in a simple manner. 
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In this case, the quantity of new air taken in the intal^e passage 
(hereinafter may be referred to as "Intake new air quantity") is 
advantageously used as the quantity of the exhaust gas discharged from the 
exhaust passage to the outside. In general, a device for measuring the 
intake new air quantity, such as airfiow meter, is provided In the intake 
passage of the engine In order to measure the quantity of new air taken in 
the intake passage. Further, in an ordinary operating state of the engine 
(in particular, in a steady operation state), the quantity of the exhaust gas 
discharged from the exhaust passage to the outside becomes substantially 
equal to the intake new air quantity. 

Accordingly, with this scheme, the quantity of NOx discharged from 
the exhaust passage to the outside can be estimated without provision of an 
additional device for measuring the quantity of exhaust gas discharged from 
the exhaust passage to the outside. As a result the NOx discharge 
quantity can be accurately estimated by use of a simple configuration. 

In the NOx discharge quantity estimation method according to the 
present invention, the combustion region Is preferably estimated by the 
steps of estimating a quantity of oxygen taken in the combustion chamber 
and a quantity of oxygen consumed as a result of combustion; and 
estimating the combustion region by use of a ratio of the quantity of oxygen 
consumed as a result of combustion to the quantity of oxygen taken in the 
combustion chamber. 

The quantity of oxygen taken in the combustion chamber can be 
obtained by multiplying the oxygen concentration (intake-gas oxygen 
concentration) of gas taken in the combustion chamber via the intake valve 
(hereinafter may be referred to as "intake gas"), which is measured by 
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means of an intake-gas oxygen concentration sensor, by tlie total quantity of 
gas tal<en in tlie combustion cliamber (lierelnafter may be referred to as 
"cylinder interior total gas quantity"). The quantity of oxygen consumed as 
a result of combustion can be determined under the assumption that all 
injected fuel completely bums completely at the stoichiometric air-fuel ratio. 

When the gas components (Including oxygen molecules and NOx) of 
intalte gas taken in the combustion chamber are assumed to be uniformly 
distributed throughout the entire region of the combustion chamber, the ratio 
of the "quantity of oxygen consumed as a result of combustion" to the 
"quantity of oxygen taken in the combustion chamber" (hereinafter referred 
to as an "oxygen quantity ratio") accurately represents the ratio of the 
volume of the combustion region to that of the combustion chamber, and 
accordingly, represents the ratio of the quantity of NOx present in the 
combustion region before combustion to the total quantity of NOx taken in 
the combustion chamber (circulated into the combustion chamber via the 
EGR apparatus). 

Accordingly, when the method is configured to estimate the 
combustion chamber by use of the oxygen quantity ratio, the 
combustion-region circulated NOx quantity (the quantity of a portion of NOx 
circulated into the combustion chamber, the portion being present within the 
combustion region before combustion) and the non-combustion-reglon 
circulated NOx quantity (the quantity of a portion of NOx circulated into the 
combustion chamber, the portion being present within the non-combustion 
region before combustion) can be obtained more accurately. 

Moreover, In the NOx discharge quantity estimation method 
according to the present invention, when the internal combustion engine is 
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configured to effect, in each operation cycle, at least one pilot injection and 
then main Injection, the combustion-generated NOx quantity (the quantity of 
NOx generated upon completion of combustion of fuel injected through the 
main injection) is preferably estimated in consideration of influences of an 
inert gas generated as a result of the pilot injection. 

In the case where main injection is performed after completion of 
pilot injection one or a plurality of times, an inert gas generated as a result 
of combustion of fuel injected through the pilot injection remains in the 
combustion region at the time when fuel injected through the main injection 
starts combustion in the combustion region, whereby the oxygen 
concentration in the combustion region decreases. Further, the highest 
combustion temperature of fuel injected through the main injection drops as 
a result of the thermal capacity of the inert gas remaining in the combustion 
chamber. 

As Is understood from above, the combustion-generated NOx 
quantity; i.e., quantity of NOx generated in the combustion region as a result 
of combustion, decreases in the case where a required amount of fuel is 
injected portionwise through pilot injection and main injection, as compared 
with the case where the required amount of fuel is injected at a single time. 

Accordingly, in the case where pilot injection is performed at least 
one time before main Injection, the above-mentioned configuration adapted 
to estimate the combustion-generated NOx quantity in consideration of 
influences of an inert gas generated as a result of the pilot injection can 
estimate the combustion-generated NOx quantity more accurately. The 
drop In the combustion-generated NOx quantity caused by the influence of 
the inert gas can be calculated as a function of, for example, fuel injection 
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quantity at the time of pilot injection or timing of the pilot injection. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 a schematic diagram showing the overall configuration of a 
system in which an engine control apparatus, which performs an NOx 
discharge quantity estimation method for an internal combustion engine 
according to an embodiment of the present invention, is applied to a 
four-cylinder internal combustion engine (diesel engine). 

FIG. 2 is a diagram schematically showing a state In which gas Is 
tal<en from an intake manifold to a certain cylinder and is then discharged to 
an exhaust manifold. 

FIG. 3 Is a flowchart showing a routine which the CPU shown in FIG. 
1 executes so as to control fuel injection quantity, etc. 

FIG. 4 is a table for determining a fuel injection quantity, to which the 
CPU shown In FIG. 1 refers during execution of the routine shown In FIG. 3. 

FIG. 5 is a table for determining a base fuel injection timing, to which 
the CPU shown in FIG. 1 refers during execution of the routine shown in FIG. 
3. 

FIG. 6 is a table for determining a base fuel injection pressure, to 
which the CPU shown in FIG. 1 refers during execution of the routine shown 
in FIG. 3. 

FIG. 7 is a table for determining a target NOx discharge quantity, to 
which the CPU shown in FIG. 1 refers during execution of the routine shown 
in FIG. 3. 

FIG. 8 is a table for determining an injection-timing correction value, 
to which the CPU shown in FIG. 1 refers during execution of the routine 
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Shown in FIG, 3, 

FIG. 9 is a flowcliart showing a routine which the CPU shown In FIG. 
1 executes so as to compute an NOx discharge quantity. 

BEST MODE FOR CARRYING OUT THE INVENTION 

With reference to the drawings, there will now be described an 
control apparatus of an Internal combustion engine (diesel engine), which 
apparatus performs an NOx discharge quantity estimation method according 
to an embodiment of the present invention. 

FIG. 1 schematically shows the entire configuration of a system in 
which such an engine control apparatus is applied to a four-cylinder internal 
combustion engine (diesel engine) 10. This system comprises an engine 
main body 20 including a fuel supply system; an intake system 30 for 
introducing gas to combustion chambers (cylinder interiors) of individual 
cylinders of the engine main body 20; an exhaust system 40 for discharging 
exhaust gas from the engine main body 20; an EGR apparatus 50 for 
performing exhaust circulation; and an electronic control apparatus 60. 

Fuel iniection valves (injection valves, injectors) 21 are disposed 
above the Individual cylinders of the engine main body 20. The fuel 
injection valves 21 are connected via a fuel line 23 to a fuel injection pump 
22 connected to an unillustrated fuel tank. The fuel injection pump 22 is 
electrically connected to the electronic control apparatus 60. In 
accordance with a drive signal from the electronic control apparatus 60 (an 
instruction signal corresponding to a base fuel injection pressure Pcrbase to 
be described later), the fuel injection pump 22 pressurizes fuel in such a 
manner that the actual injection pressure (discharge pressure) of fuel 
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becomes equal to the base fuel injection pressure Pcrbase. 

Thus, fuel pressurized to the base fuel injection pressure Pcrbase is 
supplied from the fuel injection pump 22 to the fuel Injection valves 21 . 
Moreover, the fuel injection valves 21 are electrically connected to the 
electronic control apparatus 60. In accordance with a drive signal (an 
instruction signal corresponding to a fuel injection quantity qfin to be 
described later) from the electronic control apparatus 60, each of the fuel 
injection valves 21 opens for a predetermined period of time so as to inject, 
directly to the combustion chamber of the corresponding cylinder, the fuel 
pressurized to the base fuel Injection pressure Pcrbase, in the fuel injection 
quantity qfin. 

The Intake system 30 includes an intake manifold 31, which is 
connected to the respective combustion chambers of the Individual cylinders 
of the engine main body 20; an intake pipe 32, which Is connected to an 
upstream-side branching portion of the intake manifold 31 and constitutes 
an Intake passage In cooperation with the intake manifold 31 ; a throttle 
valve 33, which Is rotatably held within the intake pipe 32; a throttle valve 
actuator 33a for rotating the throttle valve 33 in accordance with a drive 
signal from the electronic control apparatus 60; an Intercooler 34, which Is 
interposed In the Intake pipe 32 to be located on the upstream side of the 
throttle valve 33; a compressor 35a of a turbocharger 35, which is 
Interposed In the intake pipe 32 to be located on the upstream side of the 
intercooler 34; and an air cleaner 36, which is disposed at a distal end 
portion of the Intake pipe 32. 

The exhaust system 40 includes an exhaust manifold 41, which is 
connected to the individual cylinders of the engine main body 20; an 
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exhaust pipe 42, which is connected to a downstream-side merging portion 
of the exhaust manifold 41 ; a turbine 35b of the turbocharger 35 interposed 
in the exhaust pipe 42; and a diesel particulate filter (hereinafter referred to 
as "DPNR") 43, which is interposed in the exhaust pipe 42. The exhaust 
manifold 41 and the exhaust pipe 42 constitute an exhaust passage. 

The DPNR 43 is a filter unit which accommodates a filter 43a formed 
of a porous material such as cordierite and which collects, by means of a 
porous surface, the particulate matter contained In exhaust gas passing 
through the filter. In the DPNR 43, at least one metal element selected 
from ailcaline metals such as potassium K, sodium Na, lithium Li, and 
cesium Cs; aikallne-earth metals such as barium Ba and calcium Ca; and 
rare-earth metals such as lanthanum La and yttrium Y Is carried, together 
with platinum, on alumina serving as a carrier. Thus, the DPNR 43 also 
serves as a storage-reduction-type NOx catalyst unit which, after absorption 
of NOx, releases the absorbed NOx and reduces it. 

The EGR apparatus 50 includes an exhaust circulation pipe 51 , 
which forms a passage (EGR passage) for circulation of exhaust gas; an 
EGR control valve 52, which is Interposed in the exhaust circulation pipe 51 ; 
and an EGR cooler 53. The exhaust circulation pipe 51 establishes 
communication between an exhaust passage (the exhaust manifold 41) 
located on the upstream side of the turbine 35b, and an intake passage (the 
intake manifold 31) located on the downstream side of the throttle valve 33. 
The EGR control valve 52 responds to a drive signal from the electronic 
control apparatus 60 so as to change the quantity of exhaust gas to be 
circulated (exhaust-gas circulation quantity, EGR-gas flow rate). 

The electronic control apparatus 60 is a microcomputer which 
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includes a CPU 61, ROM 62, RAM 63, backup RAM 64, an interface 65, etc., 
wliich are connected to one anotiier by means of a bus. The ROM 62 
stores a program to be executed by the CPU 61 , tables (lookup tables, 
maps), constants, etc. The RAM 63 allows the CPU 61 to temporarily store 
data when necessary. The backup RAM 64 stores data In a state in which 
the power supply is on, and holds the stored data even after the power 
supply is shut off. The interface 65 contains A/D converters. 

The interface 65 is connected to a hot-wire-type airflow meter 71 , 
which serves as air flow rate (new air flow rate) measurement means, and is 
disposed In the intake pipe 32; an intake gas temperature sensor 72, which 
is provided in the intake passage to be located downstream of the throttle 
valve 33 and downstream of a point where the exhaust circulation pipe 51 is 
connected to the intake passage; an Intake pipe pressure sensor 73, which 
Is provided in the intake passage to be located downstream of the throttle 
valve 33 and downstream of the point where the exhaust circulation pipe 51 
is connected to the intake passage; a crank position sensor 74; an 
accelerator opening sensor 75; £md an intake-gas oxygen concentration 
sensor 76 provided in the intake passage to be located downstream of the 
throttle valve 33 and downstream of the point where the exhaust circulation 
pipe 51 is connected to the intake passage. The interface 65 receives 
respective signals from these sensors, and supplies the received signals to 
the CPU 61 . Further, the interface 65 Is connected to the fuel injection 
valves 21 , the fuel injection pump 22, the throttle valve actuator 33a, and 
the EGR control valve 52; and outputs corresponding drive signals to these 
components in accordance with instructions from the CPU 61 . 

The hot-wire-type airflow meter 71 measures the mass flow rate of 
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intake air (new air) passing through the intal<e passage (intake new air 
quantity per unit time), and generates a signal indicating the mass flow rate 
Ga (intake new air flow rate Ga). The intake gas temperature sensor 72 
detects the temperature of the above-mentioned intake gas, and generates 
a signal representing the intake gas temperature Tb. The intake pipe 
pressure sensor 73 measures the pressure of intake gas (i.e., intake pipe 
pressure), and generates a signal representing the intake pipe pressure Pb. 

The crank position sensor 74 detects the absolute crank angle of 
each cylinder, and generates a signal representing the crank angle CA and 
engine speed NE; i.e., rotational speed of the engine 10. The accelerator 
opening sensor 75 detects an amount by which an accelerator pedal AP is 
operated, and generates a signal representing the accelerator pedal 
operated amount Accp. The intake-gas oxygen concentration sensor 76 
detects the oxygen concentration of intake gas (i.e., intake-gas oxygen 
concentration), and a signal representing intake-gas oxygen concentration 
R02_in. 

Outline of NQv Discharge Quantity Estimation Method 

Next, there will be described an outline of an NOx discharge quantity 
estimation method according to the embodiment of the present invention 
performed by the control apparatus of the internal combustion engine having 
the above-described configuration (hereinafter may be referred to as the 
"present apparatus"). FIG. 2 is a diagram schematically showing a state in 
which gas (intake gas) is taken from the intake manifold 31 into a certain 
cylinder (cylinder interior) of the engine 10 and is then discharged to the 
exhaust manifold 41 . 

As shown in FIG. 2, intake gas (accordingly, cylinder interior gas) 
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includes new air tal<en from the tip end of the intal<e pipe 32 via the throttle 
valve 33, and EGR gas (Including NOx) taken from the exhaust circulation 
pipe 51 via the EGR control valve 52. The mass ratio (i.e., EGR ratio) of 
the mass of the taken EGR gas (EGR gas mass) to the sum of the mass of 
the taken new air (new air mass) and the mass of the taken EGR gas (EGR 
gas mass) changes depending on the opening of the throttle valve 33 and 
the opening of the EGR control valve 52, which are properly controlled by 
the electronic control apparatus 60 (CPU 61) In accordance with the 
operating condition. 

During an intake stroke, the intake gas (i.e., gas composed of the 
new air and the EGR gas containing NOx) is taken in the cylinder via an 
opened intake valve Vin as the piston moves downward, and the 
thus-produced gas mixture serves as cylinder Interior gas. The cylinder 
Interior gas Is confined within the cylinder when the intake valve Vin closes 
upon the piston having reached bottom dead center (hereinafter referred to 
as "ATDC-ISO""), and then compressed in a subsequent compression 
stroke as the piston moves upward. When the piston reaches top dead 
center (specifically, when a final fuel injection timing finjfin to be described 
later comes), the present apparatus opens the corresponding fuel injection 
valve 21 for a predetermined period of time corresponding to the fuel 
injection quantity qfin. to thereby inject fuel directly into the cylinder. As a 
result, the Injected fuel disperses in the cylinder with elapse of time, while 
mixing with the cylinder Interior gas to produce a gas mixture. The gas 
mixture starts combustion by means of self ignition at a predetermined 
timing. 

In the present embodiment, such combustion is assumed to occur 
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only In a combustion region (hereinafter may be referred to as "region B"; 
see FIG. 2), which is a portion of the combustion chamber and is estimated 
as described later, and not to occur in a non-combustion region (hereinafter 
may be referred to as "region A"; see FIG. 2), which is the remaining portion 
of the combustion chamber other than the region B. Cylinder interior gas 
remaining in the combustion chamber after combustion is discharged, as 
exhaust gas, to the exhaust manifold 41 via the exhaust valve Vout, which is 
held open during the exhaust strol<e, as the piston moves upward. The 
exhaust gas is then discharged to the outside via the exhaust pipe 42. 

Next, there will be describe a specific NOx discharge quantity 
estimation method, which is applied to the internal combustion engine 10 
equipped with the EGR apparatus 50, and which is performed by the 
present apparatus, in the NOx discharge quantity estimation method, upon 
arrival of each time when the final fuel Injection timing finjfin for a cylinder to 
which fuel is Injected (hereinafter referred to as "fuel Injection cylinder"), the 
mass of NOx contained in exhaust gas (I.e., NOx discharge quantity, actual 
NOx discharge quantity NOxact) is estimated, the exhaust gas being 
discharged from the fuel Injection cylinder to the outside via the exhaust 
valve Vout and the exhaust passage during the exhaust strol<e immediately 
after the final fuel injection timing. 

In this method, for estimation of the actual NOx discharge quantity 
NOxact, the ratio of the mass of NOx present in the above-described region 
B before combustion to the total mass of NOx taken in the combustion 
chamber (hereinafter referred to as "NOx quantity ratio RatioNOx"), the 
mass of NOx remaining In the above-described region A after the 
combustion, and the mass of NOx remaining In the above-described region 
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B after the combustion must be estimated. Tiierefore, tlie methods for 
obtaining these values will be described with reference to FIG. 2. 
<l\/!ethod of Obtaining NOx quantity Ratio RatioNOx> 

Each of gas components, including oxygen molecules and NOx, of 
intake gas taken in the combustion chamber (cylinder interior gas) is 
assumed to be uniformly distributed over the entire region within the 
combustion chamber. Further, in this state, all oxygen present within the 
region B is assumed to be consumed by combustion. In this case, the ratio 
of the "mass of oxygen consumed by combustion" to the "total mass of 
oxygen taken in the combustion chamber" (i.e., the above-mentioned 
oxygen quantity ratio) represents the ratio of the volume of the region B to 
the volume of the combustion chamber, and also represents the ratio of the 
mass of NOx present In the region B before combustion to the total mass of 
NOx taken in the combustion chamber (accordingly, the above-mentioned 
NOx quantity ratio RatioNOx). In other words, the region B can be 
estimated by use of the oxygen quantity ratio. 

Accordingly, the oxygen quantity ratio is obtained in order to obtain 
the NOx quantity ratio RatioNOx, and the "total mass of oxygen taken in the 
combustion chamber" and the "mass of oxygen consumed by combustion" 
must be obtained in order to obtain the oxygen quantity ratio. 

The "total mass of oxygen taken in the combustion chamber" can be 
obtained through muKipiication of the total mass of gas taken in the 
combustion chamber (hereinafter referred to as "cylinder interior total gas 
quantity Gcyl") by the oxygen concentration of the cylinder interior gas 
before combustion. The cylinder interior total gas quantity Gcyl can be 
obtained in accordance with Eq. (1), which is based on the state equation of 
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gasat ATDC-iaO". 

Gcyl = (PaO-VaO)/(RTaO) (1) 

In Eq. (1), PaO is bottom-dead-center cylinder interior gas pressure; 
i.e., cylinder Interior gas pressure at ATDC-1 80°. At ATDC-1 80", tlie 
cylinder interior gas pressure is considered to be substantially equal to the 
intal<e pipe pressure Pb. Therefore, the bottom-dead-center cylinder 
interior gas pressure PaO can be obtained from the Intake pipe pressure Pb 
detected by means of the intake pipe pressure sensor 73 at ATDC-1 80*. 
VaO is bottom-dead-center combustion chamber volume; i.e., combustion 
chamber volume at ATDC-1 80**. The combustion chamber volume Va can 
be represented as a function of the crank angle CA on the basis of the 
design specifications of the engine 10. Therefore, the bottom-dead-center 
combustion chamber volume VaO can be obtained on the basis of the 
function. TaO is bottom-dead-center cylinder interior gas temperature; i.e., 
cylinder interior gas temperature at ATDC-1 80°. At ATDC-1 80°, the 
cylinder Interior gas temperature is considered to be substantially equal to 
the intake gas temperature Tb. Therefore, ttie bottom-dead-center cylinder 
interior gas temperature TaO can be obtained from the intake gas 
temperature Tb detected by means of the intake gas temperature sensor 72 
at ATDC-1 80°. R Is the gas constant of the cylinder interior gas. 

Moreover, the oxygen concentration of the cylinder interior gas 
before combustion is considered to be substantially equal to the intake-gas 
oxygen concentration R02_in at ATDC-1 80°. Therefore, the oxygen 
concentration of the cylinder interior gas before combustion can be obtained 
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from the bottom-dead-center intake-gas oxygen concentration RO20Jn 
detected by means of tlie intake-gas oxygen concentration sensor 76 at 
ATDC-1 80°. From tlie above, tlie "total mass of oxygen taken in the 
combustion chamber" can be represented as "Gcyl'RO20_in." 

Meanwhile, the "mass of oxygen consumed by combustion" can be 
represented as "K Qf uel" under the assumption that the entirety of injected 
fuel (i.e., fuel in the above-mentioned fuel injection quantity qfin) burns 
completely at the stoichiometric air-fuel ratio stoich. K is a coefficient, 
which is a value obtained through multiplication of the mass ratio (0.23) of 
oxygen contained in the atmosphere by the stoichiometric air-fuel ratio 
stoich (e.g., 14.6); i.e., "0.23-stoich." Qfuel is a value equal to the 
above-mentioned fuel injection quantity qfln. From the above, the NOx 
quantity ratio RatioNOx can be obtained In accordance with the following Eq. 
(2). 



RatioNOx = (KQfuel)/(GcylRO20Jn) (2) 



<l\/lethod of Obtaining Mass of NOx Remaining in the Region A after 
Combustion> 

As described above, in the present embodiment, combustion within 
the combustion chamber is assumed to occur only in the region B and not to 
occur in the region A. Therefore, of the mass of NOx circulated into the 
combustion chamber via the EGR apparatus 50, the mass of NOx present 
within the region A before combustion (hereinafter referred to as "A-region 
circulated NOx quantity NOxA") can be considered to be conserved (held) 
within the region A as It Is, even after combustion. In other words, the 
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A-region circulated NOx quantity NOxA directly represents the "the mass of 
NOx remaining In the region A after combustion." 

The ratio of the mass of NOx present within the region A before 
combustion to the total mass of NOx taken in the combustion chamber can 
be represented by use of the above-mentioned NOx quantity ratio RatioNOx; 
I.e., represented as "1 - RatioNOx." Therefore, the A-region circulated 
NOx quantity NOxA can be obtained through multiplication of the cylinder 
interior total gas quantity Gcyl by the NOx concentration of the cylinder 
interior gas before combustion and (1 - RatioNOx). Since the NOx 
concentration of the cylinder interior gas before combustion can be 
considered to be substantially equal to the NOx concentration of intake gas 
(Intake gas NOx concentration RNOxJn), the A-region circulated NOx 
quantity NOxA can be represented by the following Eq. (3). 

NOxA = RNOxJn (1 - RatioNOx) Gcyl (3) 

In Eq. (3), the intake gas NOx concentration RNOxJn is the mass 
ratio of the mass of NOx contained in the EGR gas circulated from the EGR 
apparatus 50 to the cylinder Interior total gas quantity Gcyl. When the NOx 
concentration of EGR gas is assumed to be equal to the below-described 
exhaust gas NOx concentration RNOx_ex calculated In the previous 
operation cycle (at the time of fuel injection), the intake gas NOx 
concentration RNOxJn can be obtained in accordance with the following Eq. 
(4). 

RNOxJn = (RNOx_exGegr)/Gcyl (4) 
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In Eq. (4), Gegr is the mass of EQR gas which has been taken, as a 
portion of intake gas, from the EGR apparatus 50 into the combustion 
chamber during the intake stroke of the present operation cycle, and can be 
obtained in accordance with the following Eq. (5). 

Gegr = Gcyl - Gm (5) 

In Eq. (5), Gm represents the mass (intake new air quantity) of new 
air which has been taken, as a portion of intake gas, from the tip end of the 
intake pipe 32 into the combustion chamber during the intake stroke of the 
present operation cycle, and is calculated on the basis of the intake new air 
quantity per unit time (Intake new air flow rate Ga) measured by means of 
the airflow meter 71 , the engine speed NE based on the output of the crank 
position sensor 74, and a function f (Ga, NE) which uses the intake new air 
flow rate Ga and the engine speed NE, as arguments, so as to obtain 
quantity of intake new air per intake stroke. A bottom-dead-center intake 
new air flow rate GaO and a bottom-dead-center engine speed NEO, which 
are detected by the corresponding sensors at ATDC-180'', are used as the 
intake new air flow rate Ga and the engine speed NE, respectively. As 
described above, the A-region circulated NOx quantity NOxA, and thus the 
"mass of NOx remaining in the region A after combustion" can be obtained in 
accordance with the above-described Eq. (3). 
<Method of Obtaining Mass of NOx Remaining in the Region B after 
Combustion> 

Of the mass of NOx circulated into the combustion chamber via the 
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EGR apparatus 50, the mass of NOx present within the region B before 
combustion (hereinafter referred to as "B-region circulated NOx quantity 
NOxB2") can be represented by the following Eq. (6), as in the case of the 
A-reglon circulated NOx quantity NOxA represented by the above-described 
Eq. (3). 

NOxB2 = RNOxJn-RatioNOxGcyl (6) 

Meanwhile, as a result of combustion, NOx is generated in the region 
B (the quantity of such NOx will be referred to as "B-region 
combustion-generated NOx quantity NOxB1"). The B-region 
combustion-generated NOx quantity NOxB1 can be determined in 
accordance with the following Eq. (7), as a value obtained through 
multiplication of a combustion-generated NOx quantity per unit fuel quantity 
(hereinafter referred to as "combustion-generated NOx ratio RNOx^burn") by 
the above-mentioned fuel injection quantity Qfuel. The 
combustion-generated NOx ratio RNOx.burn is calculated on the basis of, 
for example, the fuel injection quantity Qfuel, the engine speed NE, and a 
function g(Qfuel, NE) which uses the fuel injection quantity Qfuel and the 
engine speed NE as arguments so as to obtain the combustion-generated 
NOx ratio. 

NOxB1 = RNOx_burnQfuel (7) 

As described previously, the NOx quantity (NOx concentration) of a 
gas after combustion is determined, Irrespective of NOx quantity (NOx 
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concentration) before the combustion, by means of chemical equilibrium, in 
accordance with the temperature of gas after the combustion (in such a 
manner that the higher the gas temperature, the greater the NOx quantity 
(NOx concentration)). 

Accordingly, in an ordinary operating state of the engine, in many 
cases the B-region combustion-generated NOx quantity NOxB1 (quantity of 
NOx generated within gas having an increased temperature as a result of 
combustion) becomes considerably greater than the B-region circulated NOx 
quantity NOxB2. Accordingly, basically, the B-reglon 
combustion-generated NOx quantity NOxB1 can be considered to be 
substantially equal to the "mass of NOx remaining in the region B after 
combustion"). 

However, in a case where the NOx concentration of EGR gas 
(accordingly, the below-described exhaust gas NOx concentration RNOx_ex 
calculated at the time of fuel injection in the previous operation cycle) Is 
considerably high and fuel injection quantity qfin (= Qfuel) is very small 
because of a very small accelerator operated amount Accp, the B-region 
circulated NOx quantity NOxB2 Increases, and the B-region 
combustion-generated NOx quantity NOxBI decreases because of a small 
Increase in gas temperature stemming from combustion. As a result, the 
B-reglon circulated NOx quantity NOxB2 becomes greater than the B-region 
combustion-generated NOx quantity NOxB1 in some cases. In such a case, 
the B-region circulated NOx quantity NOxB2, rather than the B-region 
combustion-generated NOx quantity NOxBI, is considered to be a more 
accurate representation of the "mass of NOx remaining in the region B after 
combustion". 
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As is understood from the above, the "mass of NOx remaining in the 
region B after combustion" can be obtained as the greater one of the 
B-region combustion-generated NOx quantity N0xB1 . which Is obtained in 
accordance with Eq. (7), and the B-region circulated NOx quantity NOxB2, 
which is obtained in accordance with Eq. (6) (the greater one will be referred 
to as "B-region final NOx quantity NOxB"). 

Once the A-region circulated NOx quantity NOxA, which is the "mass 
of NOx remaining in the region A after combustion," and the B-region final 
NOx quantity NOxB, which is the "mass of NOx remaining in the region B 
after combustion," are obtained, the total mass of NOx remaining in the 
combustion chamber after combustion can be obtained as "NOxA + NOxB." 
Moreover, the total mass of gas remaining in the combustion cheimber after 
combustion can be obtained as "Gcyl + Qfuel." 

Accordingly, the NOx concentration (exhaust gas NOx concentration 
RNOx_ex) of exhaust gas discharged from the combustion chamber to the 
exhaust passage (exhaust manifold 41) via the exhaust valve Vout during 
the exhaust stroke is equal to the mass ratio of the "total mass of NOx 
remaining in the combustion chamber after combustion" to the "total mass of 
gas remaining in the combustion chamber after combustion," and can be 
obtained In accordance with the following Eq. (8). 

RNOx_ex = (NOxA + NOxB)/(Gcyl + Qfuel) (8) 

As described above, the previous value of the exhaust gas NOx 
concentration RNOx_ex obtained in accordance with Eq. (8) Is used In the 
above-described Eq. (4) for obtaining the current value of the Intake gas 
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NOx concentration RNOxJn, under the assumption that the NOx 
concentration of exhaust gas flowing through the exhaust passage (exhaust 
manifold 41) is equal to the NOx concentration of EGR gas flowing through 
the exhaust circulation pipe 51 . 

When the NOx concentration of exhaust gas flowing through the 
exhaust passage (exhaust manifold 41 and exhaust pipe 42) Is assumed to 
be constant over the entire region of the exhaust passage, the exhaust gas 
NOx concentration RNOx_ex becomes equal to the NOx concentration of 
exhaust gas discharged to the outside from the exhaust passage 
(specifically, the end of exhaust pipe 42). 

Moreover, in an ordinary operation state (in particular, in a steady 
operation state) of the engine 10, the mass per operation cycle (exhaust 
stroke) of exhaust gas which is discharged to the outside from the exhaust 
passage (exhaust pipe 42) is substantially equal to the above-mentioned 
intake new air quantity Gm. As is understood from the above, the mass 
per operation cycle of NOx contained in exhaust gas which is discharged to 
the outside via the exhaust passage (the above-mentioned actual NOx 
discharge quantity NOxact) can be obtained in accordance with Eq. (9). 
Eq. (9) shows that as the EGR gas quantity Gegr increases, the intake new 
air quantity Gm decreases, whereby the actual NOx discharge quantity 
NOxact decreases. Accordingly, the phenomenon that the actual NOx 
discharge quantity NOxact decreases as the EGR gas quantity Gegr 
increases can be accurately expressed. 

NOxact = RNOx_exGm (9) 
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As described above, upon each arrival of tiie final fuel injection 
timing finjfin for the fuel injection cylinder, the present apparatus estimates, 
by use of Eqs. (1) to (9), the actual NOx discharge quantity NOxact; i.e., the 
mass of NOx discharged from the fuel injection cylinder via the exhaust 
valve Vout in the exhaust stroke immediately after the injection timing. The 
above is the outline of the NOx discharge quantity estimation method. 
<Outline of Fuel Injection Control> 

The present apparatus, which performs the above-mentioned NOx 
discharge quantity estimation method, calculates, at predetermined intervals, 
a target NOx discharge quantity per operation cycle NOxt on the basis of the 
above-mentioned fuel injection quantity qfin and engine speed NE. 
Subsequently, the present apparatus feedback-controls the final fuel 
injection start timing finjfin and the opening of the EGR control valve 52 in 
such a manner that the actual NOx discharge quantity NOxact estimated in 
the previous operation cycle coincides with the target NOx discharge 
quantity NOxt. 

Specifically, when the actual NOx discharge quantity NOxact 
estimated in the previous operation cycle Is greater than the target NOx 
discharge quantity NOxt, the final fuel Injection start timing finjfin to be 
applied for the fuel injection cylinder in the present operation cycle is 
delayed from the base fuel injection start timing finjbase by a predetermined 
amount, and the opening of the EGR control valve 52 is increased from the 
current degree by a predetermined amount. As a result, the highest flame 
temperature of the fuel injection cylinder in the present operation cycle is 
controlled to decrease, whereby the actual NOx discharge quantity NOxact; 
i.e., the quantity of NOx discharged from the fuel injection cylinder to the 
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outside in tlie present operation cycle, is rendered coincident witli tlie target 
NOx discharge quantity NOxt. 

Meanwhile, when the actual NOx discharge quantity NOxact 
estimated in the previous operation cycle is smaller than the target NOx 
discharge quantity NOxt, the final fuel injection start timing finjfin to be 
applied for the fuel injection cylinder In the present operation cycle is 
advanced from the base fuel injection start timing finjbase by a 
predetermined amount, and the opening of the EGR control valve 52 Is 
decreased from the current degree by a predetermined amount. As a 
result, the highest flame temperature of the fuel injection cylinder in the 
present operation cycle is controlled to increase, whereby the actual NOx 
discharge quantity NOxact; i.e., the quantity of NOx discharged from the fuel 
injection cylinder to the outside in the present operation cycle, is rendered 
coincident with the target NOx discharge quantity NOxt. The above is the 
outline of fuel injection control. 
Actual Ooeration 

Next, actual operations of the control apparatus of the internal 
combustion engine having the above-described configuration will be 
described. 

<Control of Fuel Injection Quantity, Eto 

The CPU 61 repeatedly executes, at predetermined intervals, a 
routine shown by the flowchart of FIG. 3 and adapted to control fuel injection 
quantity, etc. Therefore, when a predetermined timing has been reached, 
the CPU 61 starts the processing from step 300, and then proceeds to step 
305 so as to obtain a fuel Injection quantity qfin from an accelerator opening 
Accp, an engine speed NE, and a table (map) Mapqfin shown in FIG. 4. 
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The table Mapqfin defines the relation between accelerator opening Accp 
and engine speed NE, and fuel Injection quantity qfin; and Is stored In the 
ROM 62. 

Subsequently, the CPU 61 proceeds to step 310 so as to determine 
a base fuel Injection timing finjbase from the fuel injection quantity qfin, the 
engine speed NE, and a table Mapflnjbase shown In FIG. 5. The table 
Mapflnjbase defines the relation between fuel injection quantity qfin and 
engine speed NE, and base fuel injection timing finjbase; and Is stored In 
the ROM 62. 

Subsequently, the CPU 61 proceeds to step 31 5 so as to determine 
a base fuel Injection pressure Pcrbase from the fuel injection quantity qfin, 
the engine speed NE, and a table MapPcrbase shown In FIG. 6. The table 
MapPcrbase defines the relation between fuel injection quantity qfin and 
engine speed NE, and base fuel injection pressure Pcrbase; and is stored in 
the ROM 62. 

Subsequently, the CPU 61 proceeds to step 320 so as to determine 
a target NOx discharge quantity NOxt from the fuel Injection quantity qfin, 
the engine speed NE, and a table MapNOxt shown in FIG. 7. The table 
MapNOxt defines the relation between fuel Injection quantity qfin and engine 
speed NE, and target NOx discharge quantity NOxt; and Is stored in the 
ROM 62. 

Subsequently, the CPU 61 proceeds to step 325 so as to store, as 
an NOx discharge quantity deviation ANOx, a value obtained through 
subtraction, from the target NOx discharge quantity NOxt, of the latest actual 
NOx discharge quantity NOxact, which is computed at a fuel Injection timing 
in a previous operation cycle by a routine to be described later. 
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Subsequently, the CPU 61 proceeds to step 330 so as to determine 
an injection-timing correction value A6 from the NOx discharge quantity 
deviation ANOx and a table MapAO shown in FIG. 8. The table MapAO 
defines the relation between NOx discharge quantity deviation ANOx and 
injection-timing correction value A6, and Is stored In the ROM 62. 

Next, the CPU 61 proceeds to step 335 so as to correct the base 
fuel injection timing finjbase by the injection-timing correction value A8 to 
thereby obtain a final fuel injection timing finjfin. Thus, the fuel injection 
timing Is corrected in accordance with the NOx discharge quantity deviation 
ANOx. As is apparent from FIG. 8, when the NOx discharge quantity 
deviation ANOx is positive, the injection-timing correction value Ae becomes 
positive, and its magnitude increases with the magnitude of the NOx 
discharge quantity deviation ANOx, whereby the final fuel injection timing 
finjfin is shifted toward the advance side. When the NOx discharge quantity 
deviation ANOx is negative, the injection-timing correction value AG 
becomes negative, and its magnitude increases with the magnitude of the 
NOx discharge quantity deviation ANOx, whereby the final fuel injection 
timing finjfin is shifted toward the retard side. 

Subsequently, the CPU 61 proceeds to step 340 so as to determine 
whether the injection start timing (i.e., the final fuel injection timing finjfin) is 
reached for the fuel injection cylinder. When the CPU 61 makes a "No" 
determination in step 340, the CPU 61 proceeds directly to step 395 so as to 
end the current execution of the present routine. 

In contrast, when the CPU 61 makes a "Yes" determination in step 
340, the CPU 61 proceeds to step 345 so as to inject fuel in an amount of 
the fuel injection quantity qfin into the fuel injection cylinder from the fuel 
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injection valve 21 at tlie base fuel injection pressure Pcrbase. In the 
subsequent step 350, the CPU 61 determines whether the NOx discharge 
quantity deviation ANOx is positive. When the CPU 61 makes a "Yes" 
determination in step 350, the CPU 61 proceeds to step 355 so as to reduce 
the opening of the EGR control valve 52 from the current degree by a 
predetermined amount. Subsequently, the CPU 61 proceeds to step 370. 

When the CPU 61 makes a "No" determination in step 350, the CPU 
61 proceeds to step 360 so as to determine whether the NOx discharge 
quantity deviation ANOx Is negative. When the CPU 61 makes a "Yes" 
determination in step 360, the CPU 61 proceeds to step 365 so as to 
increase the opening of the EGR control valve 52 from the current degree by 
a predetermined amount. Subsequently, the CPU 61 proceeds to step 370. 
When the CPU 61 makes a "No" determination in step 360 (i.e., when the 
NOx discharge quantity deviation ANOx is zero), the CPU 61 proceeds to 
step 370 without changing the opening of the EGR control valve 52. 

in this manner, the opening of the EGR control valve 52 is changed 
according to the NOx discharge quantity deviation ANOx. In step 370, the 
CPU 61 stores, as the value Qfuel, the fuel injection quantity qfin actually 
injected. Subsequently, the CPU 61 proceeds to step 395 so as to end the 
current execution of the present routine. Through the above-described 
processing, control of fuel injection quantity, fuel injection timing, fuel 
Injection pressure, and opening of the EGR control valve 52 is achieved. 
<Calculation of NOx Discharge Quantjty> 

Meanwhile, the CPU 61 repeatedly executes, at predetermined 
intervals, a routine shown by the flowcharts of FIG. 9 and adapted to 
calculate actual NOx discharge quantity NOxact. Therefore, when a 
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predetermined timing has been reached, the CPU 61 starts the processing 
from step 900, and then proceeds to step 905 so as to determine whether 
the crank angle CA at the present point In time coincides with ATDC-ISO"*. 

Description will be continued under the assumption that the crank 
angle CA at the present point In time has not yet reached ATDC-180**. In 
this case, the CPU 61 makes a "No" determination in step 905, and then 
proceeds directly to step 935 so as to determine whether the fuel Injection 
start timing (I.e., the final fuel injection timing finjfin) for the fuel Injection 
cylinder has come. Since the crank angle CA at the present point in time 
has not yet reached ATDC-180'', the CPU 61 makes a "No" determination in 
step 935, and then proceeds directly to step 995 so as to end the current 
execution of the present routine. 

After that, the CPU 61 repeatedly performs the processing of steps 
900, 905, 935, and 995 until the crank angle CA reaches ATDC-180''. 
When the crank angle CA has reached ATDC-180'', the CPU 61 makes a 
"Yes" determination when it proceeds to step 905, and then proceeds to 
step 910. In step 910, the CPU 61 stores, as bottom-dead-center cylinder 
interior gas temperature TaO, bottom-dead-center cylinder interior gas 
pressure PaO, bottom-dead-center intake new air flow rate GaO, and 
bottom-dead-center engine speed NEO, respectively, the Intake gas 
temperature Tb, the intake pipe pressure Pb, the intake new air flow rate Ga, 
and the engine speed NE, which are detected by means of the intake gas 
temperature sensor 72, the intake pipe pressure sensor 73, the airflow 
meter 71 , and the crank position sensor 74, respectively, at the present 
point In time (ATDC-180^). 

Subsequently, the CPU 61 proceeds to step 915 so as to store, as 
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bottom-dead-center intake-gas oxygen concentration RO20Jn, the 
intake-gas oxygen concentration R02Jn detected by means of the 
intake-gas oxygen concentration sensor 76 at the present point in time 
(ATDC-1 80°). In the subsequent step 920, the CPU 61 computes the 
cylinder interior total gas quantity Gcyl in accordance with the 
above-described Eq. (1). Here, the values stored at step 910 are 
employed as the bottom-dead-center cylinder interior gas pressure PaO and 
the bottom-dead-center cylinder interior gas temperature TaO. 

Subsequently, the CPU 61 proceeds to step 925 so as to compute 
an intake new air quantity Gm from the bottom-dead-center intake new air 
flow rate GaO and the bottom-dead-center engine speed NEO in accordance 
with the above-defined function f. In the subsequent step 930, the CPU 61 
computes an EGR gas quantity Gegr on the basis of the cylinder interior 
total gas quantity Gcyl computed in step 920 and the intake new air quantity 
Gm, and in accordance with the above-described Eq. (5). Subsequently, 
the CPU 61 proceeds to step 935 so as to make a "No" determination, and 
then proceeds to step 995 so as to end the current execution of the present 
routine. 

After that, the CPU 61 repeatedly performs the processing of steps 
900, 905, 935, and 995 until the fuel injection timing (I.e.. the final fuel 
Injection timing finjfin) comes. When the final fuel injection timing finjfin 
comes, the CPU 61 makes a "Yes" determination in step 935 and then 
proceeds to step 940 so as to calculate the intake gas NOx concentration 
RNOxJn in accordance with the above-described Eq. (4). Here, the values 
computed in steps 930 and 920 are employed as the EGR gas quantity Gegr 
and the cylinder interior total gas quantity Gcyl, respectively. The value 
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which has been computed in step 975 (to be described later) at the fuel 
injection start timing in the previous operation cycle is employed as tiie 
exhaust gas NOx concentration RNOx_ex. 

Subsequently, the CPU 61 proceeds to step 945 so as to compute 
an A-region circulated NOx quantity NOxA in accordance with the equation 
shown in the box of step 945 of FIG. 9; the equation corresponding to the 
above-described Eqs. (2) and (3). Here, the latest value stored in step 370 
of FIG. 3 is employed as the value Qfuel. After that, the CPU 61 proceeds 
to step 950 so as to a combustion-generated NOx ratio RONx_burn on the 
basis of the value Qfuel, the engine speed NE at the present point in time, 
and the above-mentioned function g. In subsequent step 955, the CPU 61 
obtain a B-reglon combustion-generated NOx quantity NOxBI In accordance 
with the above-described Eq. (7). 

Subsequently, the CPU 61 proceeds to step 960 so as to obtain a 
B-region circulated NOx quantity N0xB2 in accordance with the equation 
described in the box of step 960 and corresponding to the above-described 
Eqs. (2) and (6). In step 965, the greater one of the B-region 
combustion-generated NOx quantity NOxBI and the B-region circulated NOx 
quantity N0xB2 is stored as a B-region final NOx quantity NOxB. 

Subsequently, the CPU 61 proceeds to step 970 so as to obtain an 
exhaust gas NOx concentration RNOx_ex in accordance with the 
above-described Eq. (8). In step 975, the CPU 61 obtains an actual NOx 
discharge quantity NOxact in accordance with the above-described Eq. (9), 
and then proceeds to step 995 so as to end the current execution of the 
present routine. After that, the CPU 61 repeatedly performs the processing 
of steps 900, 905, 935, and 995 until ATDC-180° for the fuel injection 
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cylinder comes again. 

In the above-described manner, a new actual NOx discharge 
quantity NOxact Is obtained each time the fuel injection start timing comes. 
The obtained new actual NOx discharge quantity NOxact is used in step 325 
of FIG. 3 as described above. As a result, the final fuel injection timing 
finjfin and the opening of the EGR control valve 52 to be applied to the fuel 
injection cylinder in the next operation cycle are feedback-controlled on the 
basis of the new actual NOx discharge quantity NOxact. 

As described above, in the NOx discharge quantity estimation 
method for an internal combustion engine according to the embodiment of 
the present invention, the gas components (e.g., oxygen molecules and 
NOx) of intake gas taken In the combustion chamber are assumed to be 
uniformly distributed throughout the entire region of the combustion 
chamber. Under such assumption, the combustion chamber is divided Into 
a combustion region (region B) and a non-combustion region (region A) by 
making use of the ratio of the "mass of oxygen consumed by combustion" to 
the "total mass of oxygen taken in the combustion chamber." Further, 
under the assumption that NOx generated as a result of combustion remains 
In the region B after combustion (the actual NOx quantity is the 
above-mentioned B-region final NOx quantity NOxB), and circulated NOx 
present in the region A before combustion (the quantity of circulated NOx is 
the above-mentioned A-region circulated NOx quantity NOxA) is conserved 
(remains) after combustion, the actual NOx discharge quantity NOxaxt (the 
actual quantity of NOx discharged from the exhaust passage to the outside) 
Is calculated in consideration of not only the combustion-generated NOx 
quantity in the region B but also the A-region circulated NOx quantity NOxA. 
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Accordingly, the actual NOx discharge quantity NOxaxt can be accurately 
estimated. 

The present Invention is not limited to the above-described 
embodiment, and may be modified in various manners within the scope of 
the present invention. For example, the following modifications may be 
employed. In the above-described embodiment, NOx discharge quantity 
(actual NOx discharge quantity NOxaxt) is calculated every time the fuel 
Injection start timing comes irrespective of the operating state of the engine 
10. However, the embodiment may be modified In such a manner that the 
NOx discharge quantity is calculated only when the engine is in a 
predetermined steady operation state. 

In the above-described embodiment, the above-mentioned B-region 
final NOx quantity NOxB, which is the greater one of the above-mentioned 
B-reglon combustion-generated NOx quantity NOxBI , which is obtained in 
accordance with Eq. (7), and the above-mentioned B-region circulated NOx 
quantity NOxB2, which is obtained in accordance with Eq. (6), is employed 
as the "mass of NOx remaining in the region B after combustion." However, 
the B-region combustion-generated NOx quantity N0xB1 may be always 
employed as the "mass of NOx remaining in the region B after combustion." 

In the above-described embodiment, the engine 10 is configured to 
inject fuel of a required amount (fuel Injection quantity (qfin) at a single time 
In each operation cycle. However, In the case where the engine 10 Is 
configured to inject fuel of the required amount by means of pilot injection 
and main injection in each operation cycle. In consideration of influences of 
an Inert gas generated as a result of the pilot injection, the above-mentioned 
combustion-generated NOx ratio RNOx_burn may be corrected (decreased) 
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by a predetermined amount on the basis of, for example, tlie fuel injection 
quantity at the time of the pilot injection or the timing of the pilot injection. 
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